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Abstract.  In  this  paper  we  generalize  several  results  on  uniform 
approximation  orders  with  radial  basis  functions  in  (Buhmann,  Dyn 
and  Levin,  1993)  and  (Dyn  and  Ron,  1993)  to  .^-approximation  or¬ 
ders.  These  results  apply,  in  particular,  to  approximants  from  spaces 
spanned  by  translates  of  radial  basis  functions  by  scattered  centres. 
Examples  to  which  our  results  apply  include  quasi-interpolation  and 
least-squares  approximation  from  radial  function  spaces. 


1.  Introduction 

Radial  basis  function  methods  are  tools  for  multivariable  approximation  where 
functions  /  :  R"  — »  1R  are  approximated  from  spaces  spanned  by  translates  of 
normally  just  one  function  <p :  ]Rn  — ►  ]R.  The  points  by  which  <p  is  translated 
are  usually  called  “centres”.  The  function  <p  may  or  may  not  be  a  radially 
symmetric  function,  i.e.,  p>  =  <£(||  •  Ha),  but  for  the  most  common  and  best 
studied  examples  it  is,  such  as  the  ubiquitous  multiquadric  function,  where 
<p  =  yj{-)2  +  c2,  c  being  a  parameter,  or  the  equally  well-known  thin  plate 
spline,  where  <p  =  (-)2  log(-).  This  is  why  these  methods  are  still  called  radial 
basis  function  methods,  although  their  performance  depends  much  more  on 
( p's  smoothness  and  growth  properties  than  on  its  spherical  symmetry. 

There  are  several  recent  reviews  about  the  state-of-the-art  in  the  inves¬ 
tigation  of  radial  basis  function  methods,  but  we  mention  here  only  the  one 
due  to  the  first  author  (1993).  It  is  apparent  from  this  survey  that,  in  most 
articles  about  radial  basis  functions  and  their  approximation  properties  so 
far,  only  translates  along  multiintegers  have  been  studied  in  detail.  Quasi¬ 
interpolation  and  least-squares  approximation  from  the  radial  basis  function 
spaces,  and  interpolation  at  the  centres,  have  been  investigated  thoroughly  ~ 
for  multiinteger  centres,  but,  to  the  best  of  our  knowledge,  extensions  of  this  - 
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work  to  infinitely  many  scattered  centres  have  not  been  provided  until  re¬ 
cently,  when  the  works  of  Buhmann,  Dyn  and  Levin,  and  by  Dyn  and  Ron 
(both  1993)  were  completed.  In  both  papers,  quasi-interpolation  from  radial 
basis  function  spaces  with  infinitely  many  scattered  centres  is  analyzed.  In 
the  former  paper,  quasi-interpolants  for  sets  of  scattered  centres  that  enjoy 
certain  regularity  conditions  are  explicitly  constructed.  It  is  shown  that  they 
provide  the  same  approximation  orders  as  the  quasi-interpolants  on  gridded 
centres  (which  they  agree  with  in  the  special  case  when  the  scattered  centres 
happen  to  be  gridded)  if  the  centres  become  dense  in  the  underlying  Rn.  The 
class  of  radial  basis  functions  for  which  this  approach  applies  is  restricted  by 
various  conditions,  but  none  of  the  usually  considered  fs  is  excluded. 

In  (Dyn  and  Ron,  1993),  a  different  track  is  proposed.  Rather  than 
constructing  directly  approximation  schemes  that  employ  scattered  translates 
of  if,  this  paper  suggests  and  analyzes  a  method  for  converting  known  “gridded 
centre  approximation  schemes”  to  the  scattered  case.  This  is  done  in  such 
a  way  that  the  difference  of  the  latter  and  the  former  is,  at  least,  of  the 
same  order  as  the  approximation  order  anticipated  (when  the  gridded  centres 
become  dense  in  the  underlying  space).  The  approach  to  these  results  rests 
on  the  idea  that  the  gridded  centre  approximant  can  itself  be  approximated 
by  a  function  from  the  space  spanned  by  scattered  centre  translates.  For  this, 
one  replaces  each  multiinteger  translate  of  the  radial  basis  function,  <f(-  —  j) 
for  a  j  €  ZZn  say,  by  an  approximation,  call  it  ifj,  from  the  space  spanned 
by  the  scattered  translates.  The  difference  <f(-  —j)~  ifj  then  has  to  satisfy  a 
certain  decay  condition  and  there  are  other  conditions  on  the  gridded  centre 
scheme.  This  work  applies  to  the  class  of  radial  basis  functions  as  defined 
in  (Buhmann,  Dyn  and  Levin,  1993).  It  is,  in  fact,  more  general  because 
not  just  one  quasi-interpolant  for  scattered  centres  is  manufactured.  Instead, 
a  general  conversion  method  is  established  that  allows  many  approximation 
schemes  to  be  converted  from  gridded  to  scattered  centres. 

All  the  error  estimates  in  the  literature  mentioned  so  far  are  in  T°°(Rn), 
i.e.,  they  sure  uniform  estimates.  Error  analyses  of  quasi-interpolation  schemes 
in  general  Lp(lRn)-norms  are  usually  more  subtle,  and  are  less  common  than 
their  L°°(]Rn)-counterparts.  We  mention  here  the  work  of  Lei  and  Jia  (1991), 
where  Lp(lRn  ^approximation  orders  are  established  from  spaces  spanned  by 
gridded  translates  of  a  compactly  supported  piecewise-exponential  function;  in 
this  regard  we  also  mention  arguments  employed  in  the  book  on  box  splines 
by  de  Boor,  Hollig  and  Riemenschneider  (1993)  for  a  related  theorem  on 
//^-approximation  orders.  Later,  Jia  and  Lei  (1993)  offered  approximation 
schemes  and  a  complementary  error  analysis  for  spaces  generated  by  the  grid¬ 
ded  translates  of  severed  non-compactly  supported  functions.  This  applies 
to  quasi- interpolation  with  basis  functions  of  global  support,  such  as  is  re¬ 
quired  here  for  radial  basis  function  approximants,  in  contrast  to  the  Lp- 
approximation  order  result  by  (de  Boor,  Hollig  and  Riemenschneider,  1993) 
which  only  applies  to  compactly  supported  basis  functions  if.  On  the  other 
hand,  the  scheme  proposed  by  Jia  and  Lei  is  more  specific  in  nature. 

In  our  work,  we  amalgamate  these  approaches  and  provide  Lp -error  es- 
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timates  between  the  gridded  and  the  scattered  centre  quasi-interpolants  by 
basis  functions  of  global  support.  Therefore,  they  apply  to  quasi-interpolation 
schemes  using  radial  basis  functions.  We  complement  these  results  by  sim¬ 
plified  approaches  to  estimating  the  Lp- approximation  orders  of  quasi-inter¬ 
polants  even  for  gridded  centres.  Examples  that  show  the  usefulness  of  our 
work  are  provided  in  the  paper  too. 

2.  .^-approximation  orders  on  uniform  grids 

We  always  assume  that  the  given  approximand  is  an  /  €  W£(1Rn  ),P  €  [1 ,  oo), 
without  further  mention,  IVj* (Rn)  being  the  usual  £p-Sobolev  space  of  order 
k  of  functions  in  n  variables.  We  restrict  attention  to  p  6  [l,oo)  because  the 
results  for  p  —  oo  are  readily  available,  as  pointed  out  in  the  introduction, 
and  they  are  of  the  same  nature  as  our  £p-results,  but  simpler  in  various 
respects.  In  fact,  all  our  theorems  remain  true  verbatim  if  every  occurrence  of 
p  is  replaced  by  oo.  The  positive  integer  k  is  associated  with  the  polynomial 
recovery  of  the  semi-discrete  convolution  operator 

*/>*'  :  f*-+  -j)f(j) 

j 

(all  sums  are  over  7Ln  unless  otherwise  indicated)  with  a  given  basis  function 
i.e.,  we  assume  that  x p  *'  is  the  identity  on  Iljt— l  which  is  the  space  of 
polynomials  in  n  variables  of  total  degree  less  than  k.  Such  xfis  that  are 
linear  combinations  of  radial  basis  functions  exist  in  abundance,  for  instance 
ifc  =  n  +  1  is  permissible  for  multiquadrics  and  k  =  n  +  3  for  thin  plate 
splines  when  n  is  odd  or  even  respectively  (Buhmann,  1989,  for  instance).  In 
connection  with  the  Lp-Sobolev  spaces,  we  shall  shortly  require  the  following 
further  notation.  By  |  •  |/tlP  we  shall  denote  the  homogeneous  kth  order  Lp- 
Sobolev  semi-norm  on  W£ (Rn),  i.e., 

\f\k,r  =  E 

|a|=fc 

where  D  =  (^f-,  gf~, . . . ,  ^f-)  is  the  vector  of  partial  derivatives  and  ||Dar/||J, 
is  the  Lp- norm  of  t)a  f.  Moreover,  D°  =  D°lD%2  -  •  •  D%n,  ja|  =  <*i  +  c*2  + 

•  ••  +  an  and,  as  we  shall  require  in  the  sequel,  a!  =  aila^! •  •  •  <*„!,  all  for 
a  £  2Z+  and  a  =  (ati,  a2, . . . ,  an).  Finally,  |j  •  j|P)n  and  |  •  |fc>Pin  denote  the 
£p-norm  and  the  homogeneous  kth.  order  £p-Sobolev  semi-norm,  restricted  to 
a  prescribed  subset  Q  C  Hn,  respectively,  and  ||  •  ||p,2Zn  stands  for  the  £P(Z Zn) 
norm. 

The  approximation  to  /  is  provided  in  the  form  r l>  *'  g,  for  a  suitably 
chosen  g.  In  case  /  is  smooth  enough  (say,  /  €  W* ( IR”)  H  C*(Rn)),  we  may 
take  g  =  /.  Otherwise,  a  smoothing  procedure  is  required,  and  g  is  defined 
as  g  :=  A  *  /,  where  A  €  C*(IFtn)  satisfies  the  two  decay  conditions 

HU  •  U^AHi  <  oo  and  ^  ||A||oo,>+u  <  oo. 

j 


(2.1) 
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Here,  U  :=  (— |,  j)n.  We  also  demand  A  to  satisfy  the  moment  conditions 

J  X(x)xa  dx  =  Sqc  V  |a|  <  k.  (2.2) 

Any  integrals  in  this  work  are  over  It”  unless  noted  otherwise.  We  note  that 
functions  A  that  satisfy  the  above  conditions  do  exist,  and  any  such  function 
will  do  here. 

The  approximation  scheme  /  «  0  *'  (A  *  /)  can  be  refined  by  dilation, 
that  is,  with  the  scaling  operator 

‘  f  f(-/h ), 

finer  approximations  to  /  are  provided  by 

/»**(0  *'(A  *(<Tl/hf)). 

Here  and  throughout,  *  denotes  convolution,  either  discrete  or  continuous, 
but  never  semi-discrete,  as  we  use  the  *'  notation  for  the  latter.  The  exact 
meaning  of  *  will  be  clear  from  the  context.  Note  that  this  latter  approximant 
is  a  linear  combination  of  the  /iZZn-translates  of  crh4>-  One  should  note  that 
in  order  to  derive  approximation  orders,  it  suffices  to  establish  a  bound  of  the 
form 

11/  -  0  *'  (A  *  f)\\,  <  C\f\k<p ,  (2.3) 

(where  C  denotes  here  and  everywhere  in  the  paper  a  universal  positive  con¬ 
stant),  since  then  it  is  straightforward  to  derive  the  estimate 

||/  -  <r»(*  «’  (A  .  (», /,/))) ||,  <  h  0.  (2.4) 

Our  principal  result  in  this  direction  is  as  follows: 

Theorem  1.  Let  0  :  R"  -t  R  be  a  function  that  satisfies  the  following 
conditions: 

(i)  Decay  condition:  for  all  x  6  Hn,  it  is  true  that 

W*)l  <  C(1  + 

where  k  is  a  given  positive  integer,  >  n  and  C  does  not  depend  on  x. 

(ii)  Polynomial  reproduction:  0*'  is  the  identity  on  IIfc_i. 

Let  f  €  W*(Rn),  and  let  g  €  W* (1RB)  satisfy  the  following  three  conditions: 

(a)  g  is  k-times  continuously  differentiable. 

(b)  It  is  true  that 

(E  K-j+o)  '  <  C|/k,.  (2.5) 

'  j  * 

(c)  It  is  true  that 

11/  -  *||,  <  C\f\klP.  (2.6) 
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Then 

11/  -  rfi  *'  9\\P  <  C\f\k,p.  (2.7) 

Thus,  according  to  Theorem  1,  we  choose  g  =  /,  as  long  as  /  satisfies  the 
first  two  conditions  (a)-(b)  required  of  g  in  Theorem  1.  Otherwise,  the  above 
detailed  smoothing  process  is  suitable,  i.e.,  g  :=  A  *  /  satisfies  the  three 
required  properties.  Indeed,  we  note  first  that  g  =  A* /  €  C*(Rn)fl  W*(]Rn) 
due  to  standard  properties  of  the  convolution  (see,  e.g.,  Stein  and  Weiss, 
1971),  so  property  (a)  holds  always.  The  fact  that  the  two  other  required 
properties,  i.e.,  (b)  and  (c),  also  hold  for  g  :=  A  *  /,  is  the  content  of  the 
following  two  propositions. 

Proposition  1.  Let  g  :=  A*/,  with  A  satisfying  the  moment  conditions  (2.2) 
as  well  as  the  decay  conditions  (2.1).  Then, 


11/  -  Slip  <  c l/k,, 


(2.8) 


that  is,  g  satisfies  condition  (c)  of  Theorem  1. 


Proof.  According  to  the  integral  remainder  form  of  the  multivariable  Taylor 
theorem  and  according  to  (2.2), 


(5 -/)(*)=  /*(«)  £  ^ 

J  ,n  7-  JO 


7  €Z+ 
|7|=* 


ut)dtdu ,  x6Rn. 


(2.9) 

Therefore,  using  the  generalized  Minkowski  inequality,  (2.9)  implies  that  the 
left-hand  side  of  (2.8)  is  at  most  a  constant  multiple  of 


J  jf  |A(u)||«|*(l  -  ()*-'  itdu  | /U,.  (2.10) 

By  virtue  of  (2.1),  expression  (2.10)  is  a  constant  multiple  of  \f\k,P-  □ 

The  proof  shows  that  the  constant  C  in  (2.8)  may  depend  on  A  and  k,  but  is 
independent  of  p  and  /. 

Proposition  2.  Let  g  be  as  in  Proposition  1.  Then  (2.5)  holds,  and,  further, 


\d\k,p  <  C\f\k,p- 


(2.11) 


Proof.  We  use  Holder’s  inequality  and  obtain,  for  any  f  6  Xp(lRn),  for  any 
j,  i  €  7Ln ,  and  for  every  x  €  j  +  U, 


f  *(« 

Ji+ u 


)/(i  —  u)  du 


p./-i+™l|A||,',>+u, 


(2.12) 


where  A  +  y  —  1.  Introducing  the  two  sequences 


<■(■)  l|A||p'.j+u,  HO  :=  t  €  2Z\ 
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and  summing  (2.12)  over  all  integer  t,  we  obtain  the  bound 

Therefore,  by  Young’s  inequality  (cf.,  e.g.,  Stein  and  Weiss,  1971,  p.  178), 

(£iiwii~j+uV 

'  j  ' 

Because  of  the  second  assumption  in  (2.1),  ||a||1(z»  <  oo.  Also,  ||6||Piz»  < 
2"/'||/||p,  and  we  thus  have  proved  that 

Substituting  /  :=  Dyf,  for  all  possible  7  with  I7I  =  k ,  we  obtain  the  de¬ 
sired  claim.  The  second  inequality  (2.11)  is  an  easy  consequence  of  Young’s 
inequality  and  the  fact  that,  as  a  result  of  (2.1),  A  6  Z1(Rn).  □ 

Propositions  1  and  2  assure  us  of  the  existence  of  a  suitable  g  for  the  sake  of 
constructing  the  approximation  scheme  /  «  xj>  *'  g.  We  can,  therefore,  turn 
our  attention  to  the  proof  of  Theorem  1,  as,  indeed,  we  do  now. 

The  condition  (c)  on  g  (in  Theorem  1)  allows  us  to  bound  the  error  by 

11/ -ip*1  g\\P  <  11/ - g\\P  + 1 \g-1> *' g\\P  <  C\f\k>p  + 1 \g-$  *' g\\P,  (2.13) 

and  therefore  leads  to  the  desired  estimate  (2.7),  as  soon  as  one  has  an  estimate 
of  the  form 

II?  -  ^  slip  <  CI/U.P  (2.14) 

in  hand.  One  can  in  fact  prove  more  generally  for  any  g  €  (7fc(]Rn)nkFp  (]Rn), 
and  p  €  [l,  00),  that 

\\g -</>*'  Slip  <  C\g\ktP  (2.15) 

so  long  as  g  satisfies 

ft.  (2-16) 

The  arguments  below  that  lead  to  (2.16)  also  support  the  derivation  of  (2.14) 
from  our  condition  (2.5):  one  simply  needs  to  replace  in  the  proofs  all  appli¬ 
cations  of  (2.16)  by  (2.5). 

We  first  establish  our  approximation  order  result  for  the  (relatively  simple 
case)  when  xp  is  of  compact  support.  Thereafter  we  will  provide  a  proof  for 
Theorem  1  in  its  full  generality. 
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Theorem  2.  The  conclusions  of  Theorem  1  hold  in  case  ip  is  bounded  and  of 
compact  support.  More  generally,  if  g  €  C*(R”)n  H^R”),  and  ip  is  bounded 
and  of  compact  support  and  satisfies  (ii),  then  (2.16)  implies  (2.15). 

Proof.  Let  x  €  Rn  be  arbitrary.  If  T  denotes  the  Taylor  polynomial  to  g  € 
C*(R")  about  x  of  degree  k  -  1,  we  get  (since  ip  *'  T  =  T) 


l(V»  *'$ -$)(*)!  <  ~  MM  ~  T(j)\. 


Let  m  be  a  positive  integer  that  satisfies  suppV’  C  (2m  —  1)U.  Therefore 


U*'  9-9\\i 


*[/(E 


iV’(*-;)lis(i)-^(j)l  dx 


) 


■y  i 


(2.17) 


where  Zt  :=  {j  €  TLn  \  ||r  -  j||oo  <  |(2m  -  1)}.  By  the  integral  form  of  the 
Taylor  remainder,  the  right-hand  side  of  (2.17)  is  at  most  ||V'||oo-times 


[/ (E  l»0)-T0)|)'<«* 

sj€Zx 

=  [/(E|  E  /(i  +  tj)dt 


hi-fc 


y 


dx 


7€2Z1  ^*/o  ;€Z,  J 


7€Z  + 

M=fc 


<  E  tt  /'(i-^-'llE  -<)(•)+ ty)i 


7€Zt  + 
l7l»* 


>eZ() 


where  we  have  used  the  Minkowski  inequality  and  the  generalized  Minkowski 
inequality  in  the  penultimate  and  in  the  last  inequality,  respectively. 

We  estimate  the  above  expression  term-by-term  in  7  and  for  each  t  € 
(0, 1)  separately.  Hence  we  consider 


|  £  |ZPy«l-  ()(  )  +  0)1 


(2.18) 


leaving  out  the  factor  of  (1  —  t)k  1 .  Note  that  Zz  =  Zt  if  x  €  U  +  £.  Hence 
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we  have  for  each  t  €  (0, 1)  that  (2.18)  is  at  most  the  p-th  root  of 

/ (  S  +  (1  ~  <)*)l)  dx 

=  53  /  ( 53  + (i  -  o*)i) 

= 53  /  (53 i^7^ + (!  -  *x*  -  ;))i) dx 

t  Jx]+t  \ezt  J 

-  i^lfc,oo,/+U> 


the  last  inequality  being  a  simple  consequence  of  the  fact  that  {£  +  U}/g2Zn 
form  a  tiling  of  Rn. 

Now,  if  assumption  (2.16)  on  g  holds,  we  obtain  that  the  last  expression 
is  bounded  by  a  constant  multiple  of  .  We  then  take  the  p-th  root  to 
obtain  (2.15).  If,  instead,  we  apply  condition  (b),  i.e.,  replace  the  usage  of 
(2.16)  by  using  (2.5),  we  arrive  directly  at  the  desired  result  (2.7).  □ 

We  remark  that  the  various  constants  that  appear  in  the  proof  of  Theorem  2 
depend  on  k,  m,  \\xp\\oo,  n  and  p.  These  constants,  however,  do  not  rely  on  /. 

Proof  of  Theorem  1.  Before  we  begin  the  proof  of  this  result,  we  point  out 
once  again  that  there  are  many  ips  that  satisfy  the  assumptions  of  this  theo¬ 
rem  which  are  linear  combinations  of  radial  basis  functions  (see,  for  instance, 
Buhmann,  1989,  and  see  in  particular  our  example  at  the  end  of  this  note). 

To  embark  on  the  proof,  we  first  follow  the  argument  in  the  proof  of 
Theorem  1,  until  (2.18)  is  reached.  Here,  due  to  the  global  support  of  ip,  there 
is  no  need  to  introduce  the  set  Zt.  Instead,  we  use  the  decay  assumption  (i) 
on  ip,  and  obtain  the  bound 


IIV'  *’  9~  g\\P  ^ 


w*-;)iij(;)-ro')i)P<fc] 


1 

P 


7  €7Ln+  j 

M =* 


dt, 


where  we  have  used  Minkowski’s  inequality  and  its  generalization  as  in  the 
previous  proof.  We  observe  that  the  integrand  of  this  expression  is  for  each 
f  €  ( 0,1)  less  than  or  equal  to 


-illoo) 


-Vi, 


+  (1  -  <)(■))! 


p 

p,U+i 


(2.19) 
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(ElDi+i-ii-r** 

^  «  j 

^  fe(Di +■>«•)■•* 

^  i  '  j 


\D^gii  +  tj  +  (1 


I Vg(i  +  tj  +  (1  -  t)(0)l 


||DM*-+o-  + OIL, 


(2.20) 


where  the  last  estimate  depends  on  the  triangle  inequality  for  the  uniform 
norm  and  on  t  €  (0, 1). 

For  t  =  1,  the  last  expression  can  be  identified  as  ||2?  *  G ||p;z«,  with 

B(j)  :=  (1  +  IUIU)-*,  j  £  ZZ”,  and  GO)  :=  ||JTjO  +  0|L,U.  J  €  Z!”.  For 

that  case,  Young’s  inequality  can  be  invoked  to  yield 


\\B*G\\Pt7Ln  <  \\B\\i^n\\G\\Pxn • 


While  for  t  <  1  the  expression  in  question  is  not  a  convolution  product,  it  can 
still  be  majorised  by  an  appropriate  convolution,  as  the  following  argument 
shows.  For  each  t  6  (0, 1)  fixed,  we  define  the  map  N(j)  :=  [(1  —  t)j\ ,  j  6  7L'\ 
with  [-J  denoting  (componentwise)  the  greatest  integer  less  than  (•  +  |).  We 
can  then  write  the  above  sequence  B  as  the  sum  of  sequences  of  disjoint 

supports  defined  as  follows: 


t  =  N(j), 

otherwise, 


jJe7Zn. 


By  using  the  bound 

+  tj  +  )|L,u  ^  ||GTs(i  +  j  -  N(j)  +  OIL,*,, 


defining  G(j)  :=  || D1g(j  +  OH^u’  J  €  ZZ" ’  and  using  the  nonneSativity  °f 
the  entries  of  the  sequences  B  and  G,  we  can  bound  the  last  expression  in  the 
chain  of  inequalities  (2.20)  from  above  by 


where  El  :  /  *-*■  /(•  —  ^).  Therefore,  Young’s  inequality  applies  once  again  to 
yield  the  upper  bound 


PlkznHGIUn  = 


X)U +  «»)"*(£  cor) 


JL 

p 
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Finally,  the  sum  (£;  G(j)p)l^v  ean  certainly  be  bounded  by  a  constant  multi¬ 
ple  of  (J2 j  G(j)p)1/P,  and  this  latter  expression  is  bounded  by  C|/|iiP,  by  our 
assumption  (b)  (or  by  C\g\k,p,  if,  instead  of  assuming  (b),  we  adopt  (2.16)). 
The  theorem  is  proved.  □ 

If  we  apply  Propositions  1  and  2,  we  arrive  at  the  following 

Corollary  1.  If  f  6  W*(lRn),  if  the  assumptions  (i)-(ii)  of  Theorem  1  hold 
and  if  A  €  C*(Hn)  satisfies  (2.2)  and  (2.1),  then  (2.3),  and  thus  (2.4)  too, 
hold.  Furthermore,  if  f  is  k-times  continuously  differentiable  and  satisfies 

(£m.-j+u)'sci/u,„ 

'  j  ' 

then 

»'  «ukf)h  <  Cfc*|/U,, 

for  h  tending  to  zero.  □ 

The  second  assertion  in  the  above  theorem  says  in  particular  that  for  suffi¬ 
ciently  smooth  /  one  needs  to  know  only  the  values  f\h.7in  in  order  to  compute 
its  Lp-approximant.  This  is  especially  important  from  a  practical  point  of 
view,  since  computing  the  values  on  H7Zn  of  A  *  f  can  be  hard. 

Our  final  note  in  this  section  concerns  the  conditions  imposed  on  xp  in 
Theorem  1.  The  decay  condition  (viz.,  (i))  is  essential  to  the  approximation 
orders  derived.  Though  approximation  orders  can  be  derived  under  weaker 
decay  conditions,  such  derivations  either  involve  some  complement  ary  infor¬ 
mation  of  the  exact  decay  rates  of  xp  (as  in  Buhmann,  1989,  or  in  Dyn,  Jackson, 
Levin  and  Ron,  1992),  or  require  a  substantial  modification  of  the  approxi¬ 
mation  schemes  as  well  as  their  error  analyses  (as  is  the  case  in  (de  Boor  and 
Ron,  1992a)  and  (de  Boor,  DeVore  and  Ron,  1993)).  However,  the  second 
condition,  i.e.,  that  xp*1  be  the  identity  on  II*_i  can  be  relaxed  to  (the  more 
standard  one)  that  xp*1  induces  a  linear  isomorphism  on  n*_i .  An  appropriate 
approximation  scheme  in  this  case  has  the  form 

/  «  0  *'  (p  *  A  *  /), 

with  A,  as  before,  satisfying  (2.2)  and  (2.1),  and  p,  say,  compactly  supported 
and  satisfying 

tp*’(p*p)=p,  p  €  Ilfc— i  - 

Such  constructions  are  discussed  in  detail  in  the  literature,  and  we  refer  the 
reader  to  de  Boor  and  Ron  (1992b)  where  such  schemes  were  thoroughly  anal¬ 
ysed  for  compactly  supported  xp.  The  modifications  required  in  our  arguments 
in  order  to  cover  that  more  general  setup  are  mostly  straightforward.  More¬ 
over,  one  should  notice  that,  if  p  is  selected  to  have  its  support  lying  in  7Ln , 
then  xp  *'  (p  *  •)  =  (xp  *  p )*',  and,  for  the  new  function  $  :=  xp  *  p,  is  the 
identity  on  II*_i .  Therefore,  our  results  apply  directly  to  this  important  case. 
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3.  Z^-approximation  orders  using  scattered  centres 

We  use  the  results  of  the  previous  section  in  the  derivation  of  approximation 
schemes  that  employ  scattered  translates  of  the  basis  function  <p.  The  ap¬ 
proach  as  well  as  the  error  analysis  follow  the  work  of  (Dyn  and  Ron,  1993), 
where  Z°°-approximation  orders  were  olcained. 

We  let  xl>  be  a  linear  combination 

iJ>  =  V*'n,  (3-1) 

where  is  a  given  function  in  C(lRn)  such  as  one  of  our  radial  basis  func¬ 
tions  and  where  fj,  :  7Zn  — *  1R  is  such  that  the  infinite  sum  in  (3.1)  con¬ 
verges  uniformly  on  compacta.  We  shall  as  before  consider  the  approxima¬ 
tion  operator  /  Lg  :=  rp  *'  g ,  where  g  =  X  *  f  with  A  €  Cfc(lRn)  sat¬ 
isfying  (2.2)  and  (2.1).  We  make  assumptions  according  to  those  in  (Dyn 
and  Ron,  1993):  we  always  assume  from  now  on  that  the  thus  constructed 
L  :  C(lRn)  D  X°°(lRn)  — >  C(lRn)  D  X°°(Rn)  is  a  bounded  operator  without 
further  mention.  As  in  the  previous  section,  we  use  here  stationary  schemes, 
i.e.,  define  the  fi-refinement  Lh  by  dilation: 


Lh  :=  OhLai/h. 


We  shall  further  assume  that 


KOI  <  C(1  +  IKIloo)-^,  i  6  TLn , 

W*)l<C(i  +  |MU)--,  <sk‘,  " 


(3.2) 


We  take,  for  a  given  set  of  (scattered)  centres  E  C  R",  approximations  to 
each  translate  <p(-  —  j)  of  the  form 


<pj 

«€H 


where  {Ajtz}j£2Zn,ze=.  are  suitable  coefficients  which  provide  that 

M*)  -  -j) I  <  C(  1  +  III  -  jIU)-"-  (3.3) 

holds  for  a  >  n  and  a  positive  constant  C  which  is  uniform  in  x  and  j.  We 
call  the  functions  {<fj}j  pseudo-shifts  (as  they  are  approximants  to  the  “true” 
shifts  {v?(-  —  j)}j  from  our  space),  and  use  these  pseudo-shifts  to  construct 
another  set  of  approximations,  now  to  the  localized  functions  {0(-  —  j)}j- 
We  carry  out  this  second  task  by  imitating  the  localization  process  of  (3.1). 
Precisely,  we  define 

:= 

t 


(3.4) 
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Finally,  the  approximation  scheme  LA,  that  employs  the  scattered  translates 
—  f)}*€S)  is  obtained  from  L  by  replacing  each  tp(-  —  j)  by  xpj: 

La  :f*-* 

j 

As  before,  there  are  functions  /  for  which  convolution  with  A  need  be  per¬ 
formed: 

/  «  La( A  *  /). 

The  /i-refinement  here  is  the  same  as  before  for  L,  that  is 


Lh.A  ■=  0hLAO i/fc. 

One  verifies  that  Lh,A  maps  into  the  span  of  the  /iE-translates  of  We 
remark  that  other,  more  flexible  refinements  may  be  allowed,  too,  and  refer 
the  interested  reader  to  the  discussion  of  that  point  in  §2.6  of  (Dyn  and  Ron, 
1993)  and  to  the  specific  strategy  suggested  for  the  purpose  of  treating  sets 
of  scattered  centres  becoming  dense  in  the  underlying  space  in  the  paper 
(Buhmann,  Dyn  and  Levin,  1993). 

Theorem  3.  Let  <p,  p,  xj;,  tfrj,  £*  and  Lk,a  be  as  above.  Assume  that  (3.2) 
and  (3.3)  hold.  Suppose  also  that 

Vn>n  +  k  +  e-  l  for  an  e  €  (0, 1]  ^3.5) 

and 

p  :  f  y satisfies  pUk-i  =  {0}.  (3.6) 

t 

Then,  for  h  €  (0, 1)  and  f  €  Wp (2Rn), 

ll(i*-IklM)(A./)||,<CA‘(|/|t-I,1,  +  l/k,){^A  <3-7> 

Moreover,  if  e  =  1  and  j/z(j ) | |j( <  oo,  then 

||(L*  -  L,,a)(\  *  nil,  <  Cft*|/1*,,.  (3.8) 

Implicit  in  the  statement  of  theorem  is  the  assertion  that  La{ A  *  /)  is  well- 
defined,  i.e.,  that  the  series  LA( A  *  /)  converges  absolutely  for  every,  say, 
bounded  /.  This  well-definedness  has  been  established  in  Dyn  and  Ron  (1993); 
cf.  Lemma  2.2.6  in  that  article. 

Proof.  We  begin  with  establishing  (3.8).  We  may  consider  L  and  LA  instead 
of  Lh  and  Lh,A,  because,  as  we  have  pointed  out  earlier,  the  powers  of  h  are 
introduced  by  a  suitable  scaling.  Set  g  :=  A  *  /.  Since  the  series  that  define 
L  and  La  converge  absolutely,  we  may  sum  (L  —  LA)g  by  parts  to  obtain 

\\(L-LA)g\\p  <  +^)|  • 

j  t  P 
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Since  ft  annihilates  II*_i  by  (3.6),  we  have,  for  each  j  €  ZZn, 
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E  K-t)nU  +  0  =  E  +  0  -  W  +  £)), 

e  t 

with  Tj  the  ( k  —  l)-degree  Taylor  polynomial  of  g  about  j.  Using,  as  in  the 
previous  section,  the  integral  remainder  formula  for  the  Taylor  expansion, 
followed  by  the  triangle  inequality  and  summation  by  parts,  we  obtain  the 
following  upper  bound  on  |J(X  —  L,A)g\\p : 

t  7€zz;  Ja  j  p 

hi =* 

Using  the  Minkowski  inequality  we  can  take  the  summations  over  i  and  7 
01  «  of  the  p-norm  expression.  Furthermore,  using  the  generalized  Minkowski 
inequality  we  can  also  interchange  integration  over  t  and  the  Lp- norm.  This 
leaves  us  with  the  task  of  estimating  the  expression 

||£>(-  -  J)  -  +(1  -  W|  |  •  (3.9) 

j 

Because  of  (3.3),  we  further  can  bound  expression  (3.9)  by  a  constant  multiple 

°f  „ 

||][>  +  II  •  Tilloo)“^  \D^g(j  +  (1  -  00|  (310) 

j 

The  estimation  of  expression  (3.10)  follows  a  similar  estimation  done  in  the 
proof  of  Theorem  1,  only  the  present  case  is  in  fact  slightly  simpler.  Indeed, 

|E3(h-ii  •  -jo*,)-- 1  d'9u + (i  -  wi||' 

j 

=  ElB1 + » ■  + (i  - wi||' 

»  j 

<  Etoi  +  IUIU)"'',l^(i-(  +  (i  -  wiV 

*  '  J  ' 

<cEi»’j(i +(i-wr 
1 

-  ^E  ^lit.oc.U-H 
i 

<  c|/K^ 


where  we  have  used  the  Young  inequality  in  the  second  estimate  and  Propo¬ 
sition  2  in  the  last.  Now  we  may  use  the  fact  that  Yhi  |p(— f)||f|fc  *s  finite 
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and  take  p-th  roots  to  deduce  the  estimate  (3.8)  from  the  above  chain  of 
inequalities. 

The  estimate  (3.7)  is  proved  similarly,  but  we  can  no  longer  rely  on  the 
convergence  of  ^2e  Ia*(  — £)l!  li  need  to  split  the  sum  over  £  into  two  partied 
sums:  and  the  remaining  sum.  The  first  sum  is  estimated  exactly 

as  above,  and  the  bound  obtained  then  takes  the  form 

ci/U,  £  W-<)l|£| 

t€7Znnh-'V  '' 

for  each  7  with  ('y |  =  k.  In  the  remaining  sum,  the  zero  at  j  of  g(j  +  £)  — 
T](j  +  £),  though  being  of  order  ky  is  treated  as  being  of  order  k  —  1,  thus 
leading  to  a  Taylor  remainder  formula  of  the  form 

£  ['  £  w-  -  j)  -  wi|o’(s  -  ms + a  -  «k)|  a. 

-r€2Z;  Jo  7'  j 

M=*-i 

Since  D^Tj  =  D1  g(j)  for  all  admissible  7,  we  can  follow  the  argument  that 
was  used  to  prove  (3.8),  to  obtain  a  bound  of  the  form 

£  K-<)l|£ 

t€Zn,tgh-lV  7’ 

for  each  7  with  I7I  =  k  —  1.  (In  the  derivation  of  this  latter  estimate,  one 
invokes  Proposition  2,  with  k  there  replaced  by  k  —  1.)  Now,  we  can  use  the 
fact,  that,  whenever  the  required  decay  estimate  on  p  holds  for  e  €  (0, 1), 

£  W-0l|§|  <  ch-\  (3.11) 

*€2Z"n/i-iU  7' 

£  W-«)l|£|  <<?/.*,  (3.12) 

whereas  if  e  =  1,  the  right-hand  side  of  (3.11)  has  to  be  replaced  by  — Clog/i. 
All  these  estimates  were  proved  for  instance  in  (Buhmann,  Dyn  and  Levin, 
1993),  and  are  also  well-docuntfShted  in  Dyn  and  Ron  (1993).  In  summary, 
the  total  error  bound  for,  say,  e  <  1  takes  the  form 

||(£  -  La)( A  .  /)||p  <  C(ft'-‘  I/I*,,  +  fc'|/U-i,)-  (3.13) 

It  is  straightforward  to  insert  the  log  h-tevm  into  this  expression  when  t  =  1. 
When  finally  estimating  ( Lh  —  Lh,A){ A  *  /),  we  simply  invoke  (3.13)  with  / 
replaced  by  for  the  given  value  of  h,  and  then  scale  to  obtain  (3.7).  □ 

In  analogy  to  Theorem  2.3.1  in  (Dyn  and  Ron,  1993)  (where  the  case 
p  =  00  was  discussed)  we  establish  the  following  result.  Its  purpose  is  to  relax 
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the  decay  conditions  on  p  in  the  previous  result  which  requires  p  to  decay  fast 
enough  so  as  to  be  summable  against  polynomials.  It  is  now  replaced  by  a 
condition  on  the  Fourier  series  of  p,  a  condition  that  may  be  valid  even  when 
p  decays  in  a  milder  fashion.  Incidentally,  we  will  see  that  the  rates  derived 
in  that  result  may  be  superior  to  those  of  Theorem  3. 

In  order  to  present  the  next  theorem,  we  have  to  introduce  some  addi¬ 
tional  notation.  Given  1  <  p  <  oo,  let  Wp  (Rn)  be  the  space  of  all  functions  / 

whose  Fourier  transform  /  satisfies  (1  +  |j  ■  ||2)fc/  €  Lp  (Rn),  where  j  +  jjr  =  1. 

The  Lp  -norm  of  ||  •  || */  is  denoted  by  \f\'k  .  Note  that,  for  a  nonnegative 
integer  k  and  p  >  2,  the  Hausdorff- Young  Theorem  shows  that  the  space 
W*(lRn)  equipped  with  the  norm  ||/U*  p  :=  ||(1  4- 1|  •  ||2)*/j|p'  is  continuously 
embedded  into  Wk (]Rn).  Further,  for  p  =  2  and  any  nonnegative  integer  k, 
the  two  spaces  coincide.  We  have  the  following 

Theorem  4.  Assume  1  <  p  <  oo.  Let  a 11  the  assumptions  of  the  previous 
theorem  be  in  force  (including  in  particular  (3.2))  except  for  (3.5)  and  (3.6). 
Let  f  €  Wp(lRn)  be  bounded  and  assume  that  the  inequality 

\pf\ <  C\f\'kiP  (3.14) 

holds.  Then  the  estimate 

11(1  -  La)/ II,  <  C\f\'k,p  (3.15) 

holds  as  well.  The  estimate  (3.15)  remains  valid  even  if  we  replace  f  on  its 
left-hand  side  by  (A  *  /),  provided  that ,  in  addition  to  all  other  requirements 
on  A,  A  €  L°°(Rn). 

Note  that  the  dash-semi-norm  |/|'t  reacts  to  scaling  in  the  correct  manner: 
upon  scaling  of  approximand  and  approximant  suitably  by  h,  one  does  get  in 
(3.15)  the  power  hk  into  the  bound,  as  required.  Note  also  that  the  smoothing 
process  /  «  A*  /  is  permitted  but  not  required  in  Theorem  4.  It  may,  however, 
actually  be  required  for  guaranteeing  that  the  gridded  centre  approximant  to 
the  function  /  approximates  it  to  the  best  order  established  in  Theorem  1. 

We  will  discuss  the  nature  of  condition  (3.14)  after  proving  this  theorem. 
In  particular,  we  will  compare  that  condition  to  the  polynomial  annihilation 
condition  (3.6)  assumed  in  Theorem  3. 

Proof  of  Theorem  4 •  The  proof  follows  the  argument  used  in  Theorem  3,  but 
is  much  simpler.  In  Theorem  3,  we  bounded  the  error  by 

£  !*>(■  -i)  -  wl|2>-Q/U  +  4  (3.16) 

i  t 

Since  no  smoothing  is  required,  we  replaced  g  by  /  in  (3.16).  We  bound  the 
error  in  the  same  fashion  in  the  present  proof,  and  remark  that,  since  /  is 
bounded  and  p  €  i\(7Ln)  by  (3.2),  the  sum  p(f( ■  +  j))  =  +  0 
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converges  absolutely.  By  our  assumption  (3.14)  in  the  theorem  and  the  fact 
that  the  semi-norm  |  •  |'fc  is  translation-invariant,  we  have  | p(f(-  +  j))|  < 
C\f\'k  .  We  thus  obtain  tne  bound 

j 

The  desired  result  is  now  a  consequence  of  the  assumption  (3.3).  The  proof 
so  far  also  provides  that 

||(I  -  La)( a  .  /)||,  <  C|A  .  /Ii,,  <  C||A|U|/|iiP,  (3.17) 

thus  showing  that  the  difference  (L  —  La)(^  *  f)  is  bounded  to  the  asserted 
order  too.  The  first  inequality  in  (3.17)  is  a  consequence  of  (3.15)  and  the 
second  one  is  a  consequence  of  the  definition  of  our  dash-semi-norm.  The 
theorem  is  proved.  □ 

Of  course,  Theorem  4  is  useful  only  if  one  can  establish  verifiable  condi¬ 
tions  on  /  which  ensure  (3.14).  In  this  regard,  we  remind  the  reader  that,  if 
/  fails  to  satisfy  the  conditions  of  the  theorem,  for  instance  because  /  is  not 
bounded,  we  may  still  invoke  the  result  with  respect  to  a  mollified  g  :=  A  *  /, 
and  combine  (3.15)  with  an  inequality 

Ini;.,  <  Cl/lt.,  (3.18) 

to  obtain  the  desired  estimate.  This  is  of  course  what  we  already  did  at  the 
end  of  the  above  proof,  where  (3.18)  was  a  consequence  of  the  boundedness 
of  A’s  Fourier  transform. 

The  crucial  bound  (3.14)  holds  under  a  variety  of  possible  conditions  on 
/,  k  and  /z,  and  it  is  beyond  the  scope  of  this  article  to  discuss  all  possible 
versions.  The  following  seems  to  be  the  most  useful  variant: 

Proposition  3.  Assume  /  is  bounded  and  f  6  L^R")  (~1  W* (Rn).  If 

||  •  Hi"*?  €  £'(IT),  (3.19) 

/z  being  the  Fourier  series  of  n,  then  (3.14)  holds. 

Proof.  Recall  that  /x  €  t?i(ZZn),  which  renders  \if  well-defined.  The  crux  in 
the  proof  is  the  identity 

pf  =  J  Pf-  (3.20) 

This  can  be  derived  as  in  (Dyn  and  Ron,  1993)  (cf.  the  proof  of  the  lemma  in 
Theorem  2.3.1  in  that  article),  where  the  Poisson  Summation  Formula  (Stein 
and  Weiss,  1971,  p.  252)  is  the  salient  ingredient  to  the  proof  and  which  else 
requires  only  /  €  L1(Mn)  and  boundedness  of  /.  How  it  follows  from  an 
application  of  Holder’s  inequality  tha‘ 

M  <  (/ (II  ■  K*?)') '  (/ (II  ■  115/)”’)  *  •  (3.2i) 
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Combining  (3.19)  and  (3.21)  leads  to  (3.14).  □ 

/V 

The  assumption  /  €  L1(JRn)  is  not  very  restrictive.  Indeed,  for  k  >  n/p 
that  condition  is  implied  by  the  requirement  /  €  W*(Rn).  When  k  is  too 

small  and  /  £  £1(lRn),  we  still  have  A  *  /  €  JL1^”)  if  A  is  sufficiently  smooth. 
This  smoothness  requirement  of  A  can  be  shown  to  be  mild,  and  it  decreases 
with  the  increase  of  p.  For  example,  if  p  >  2,  we  only  need  A  to  be  in  £2(Rn) 
which  we  assume  anyway.  In  short,  for  an  appropriately  chosen  A  and  under 
the  conditions  assumed  in  Theorem  4,  save  the  boundedness  of  /  and  the 
inequality  (3.14),  the  bound 

||(I  -  Ia)(A  *  /)||,  <  C\f\'kiP  (3.22) 

is  valid,  as  soon  as  (3.19)  holds. 

We  see  that  the  important  assumption  in  the  proposition  is  the  one  with 
respect  to  /i,  and  we  would  like  to  compare  it  to  the  polynomial  annihilation 
condition  (3.6)  assumed  in  Theorem  3.  If  p  decays  fast  enough  to  satisfy, 
say,  the  condition  |/i(j  )||_/ )*  <  oo  which  was  assumed  in  Theorem  3,  then 
p  is  k- times  continuously  differentiable  everywhere  and  in  particular  at  the 
origin.  It  follows  easily  that  (3.6)  implies  that  p  has  a  zero  of  order  k  at 
the  origin.  If  k  >  n/p ,  this  certainly  implies  (3.19).  Moreover,  this  latter 
condition  is  valid  for  ||  •  ||jTA  p ,  with  k'  <  k  +  n/p.  By  choosing  n/p  <  k'  < 
k  +  n/p,  we  therefore  obtain  from  Theorem  4  the  bound  (3.22),  which  yields 
that,  upon  scaling,  the  difference  between  the  scattered  centre  approximant 
and  the  gridded  approximant  decays  to  zero  with  order  hk  .  This  is  very 
satisfactory,  since  the  approximation  order  we  expect  Lh  to  provide  (in  the 
best  circumstances)  is  only  k  (not  to  mention  the  fact  that  we  had  relaxed 
the  decay  condition  on  p).  On  the  other  hand,  we  should  mention  that,  for 
example,  when  p>  2,  the  space  of  approximands  to  which  Theorem  4  applies 
(viz.,  W*(Rn))  is  only  a  proper  subspace  of  the  Sobolev  space  Wk{ Rn)  to 
which  Theorem  3  applies. 

As  an  application  of  our  results,  we  consider  fhe  orthogonal  projection 
on  a  space  spanned  by  translates  of  a  radial  basis  function.  A  comprehensive 
analysis  of  the  best  least  squares  approximations  onto  radial  functions  spaces 
and  their  approximation  orders  can  be  found  in  Ron  (1992),  but  we  will 
use  our  theorems  in  this  paper  to  deduce  £p-convergence  orders  of  best  least 
squares  approximations.  We  take  9  to  belong  to  a  class  specified  in  (Buhmann 
and  Micchelli,  1992)  which  contains  all  the  functions  that  are  mentioned  in 
the  current  literature,  see,  e.g.,  the  survey  by  Buhmann  (1993).  We  do  not 
state  here  the  full  set  of  assumptions  required  for  9  to  be  in  that  class,  but 
mention  the  salient  condition  that  9  must  have  a  singularity  at  the  origin, 
viz.,  <p(t)  ~  ||*|irfc+1-j  f°r  small  ||t||2,  where  S  €  (0,1]  and  k  is  a  positive 
integer.  It  must  also  satisfy  a  decay  condition  |9(t)|  =  0(||<||^n-e)  for  e  >  0 
and  ||f||2  — ►  00.  The  space  of  approximants  is  then  defined  as 
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S2(<f)  •-  -j)  | 

^  j 

\dj\  <  for  some  positive  £'|  D  £2(Rn). 

In  fact,  spaces  other  than  the  above  Si(<p)  can  be  used  as  the  underlying 
approximation  space  but  we  restrict  attention  to  S2(<p)  in  this  example  for 
the  sake  of  simplicity.  We  have  set  the  notations  in  such  a  way  that  k  here 
can  be  identified  with  the  k  of  Theorem  1.  The  approximation  scheme  we  use 
is  the  orthogonal  projector  P  :  L2(lRn)  — ►  S2(ip).  It  follows  from  the  required 
properties  of  ip  that  P  can  be  written  as 

PI  =  Y,UM--m(-i),  /ei2(R"),  (3.23) 

j 

as  we  shall  explain  in  the  sequel.  The  inner  product  in  (3.23)  is  the  standard 
inner  product 

j  fJtt  fit  ft  €  Z2(R").  (3.24) 

Further,  the  function  tp  in  (3.23)  is  a  function  in  S2(<f)  whose  multiinteger 
translates  are  mutually  orthonormal  with  respect  to  (3.24)  and  which  can  be 
defined  by  its  Fourier  transform 

$  =  .  g  (3.25) 

tjEj  («•  - 

The  orthonormality  property  we  claim  tp  to  have  is  easily  verified.  Further, 
it  is  established  in  Theorem  3.1  on  p.  330  and  Theorem  3.3  on  p.  333  of  the 
same  paper  by  the  first  author  and  Micchelli  that  for  <p  from  the  said  class 
of  radial  basis  functions,  tp  satisfies  the  decay  assumption  on  tp  imposed  in 
the  statement  of  Theorem  1  for  =  n  +  8.  It  therefore  decays  fast  enough 
to  admit  actually  any  /  €  i^(Rn),  p  €  [l,oo],  in  (3.23),  namely,  for  every 
/  €  Lp(Hn ),  the  series  in  (3.23)  converges  unconditionally  and  uniformly  on 
compacta  to  an  Lp( Rn)  function.  The  sum  in  (3.23)  can  be  rearranged  so 
that  the  coefficients  associated  with  the  multiinteger  translates  of  the  radial 
function  <p  decay  fast  enough  to  satisfy  the  requirement  in  the  definition  of 
our  space  S2(<p)y  cf.  Theorem  3.2  on  p.  332  in  the  mentioned  paper.  Therefore 
it  is  settled  that  (3.23)  is  indeed  the  desired  projection  operator. 

We  need  to  verify  that  the  remaining  conditions  on  tp  of  Theorem  1  and  3 
hold  and  that  (3.23)  is  actually  of  the  form  of  the  approximation  operators  we 
study  in  this  work.  To  begin  with,  we  observe  that  the  coefficients  p  that  form 
tp  as  the  semi-discrete  convolution  tp  =  <p  *'  p  satisfy  the  decay  assumptions 
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of  Theorem  3  for  e  =  6,  and  the  singularity  of  9  at  the  origin  implies  (3.6) 
because 

Next,  we  need  to  know  that  the  coefficients  of  ip(-  —  j),  j  €  2Zn,  in  the 
present  approximation  scheme  (3.23)  are  derived  from  a  convolution  process 
with  an  appropriate  convolutor  A.  Indeed,  taking  A  :=  ip( — ),  we  certainly 
have  (f,ip(-  —  j ))  =  (A  *  f)(j).  Also,  as  proved  in  Theorem  3.1  on  p.  330  of 
the  mentioned  paper,  ip  decays  in  such  a  way  to  provide  (2.1).  It  satisfies  the 
moment  conditions  (2.2)  for  the  present  k  due  to  <J?’s  singularity  at  the  origin 

/s 

and  due  to  ip’s  form  (3.25).  Consequently,  Theorems  1  and  3  are  applicable 
to  the  present  case  as  soon  as  we  show  that  ip*'  is  the  identity  on  II*_i. 
That  final  requirement  follows,  by  a  standard  argument,  from  the  fact  that  ip 
satisfies  the  so-called  Strang-Fix  conditions,  namely 

D*f{ 2*0)  =  6006*0,  p  €  ZZ",  | oc |  <  k. 

These  conditions  can  be  verified  directly  from  the  definition  of  ip  in  (3.25) 
and  from  the  singularity  of  the  radial  function’s  Fourier  transform  at  the 
origin.  Therefore,  our  Theorems  1  and  3  apply  to  this  example  and  provide 
the  desired  .^-approximation  orders. 
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